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Abstract: Mechanisms for the observed large Br{B — >■ 77' + Xs) 
are examined. We propose that the dominant fraction of the 
B f]' + Xg rate is due mainly to b — > sg*, where g* is an 
off-shell gluon, followed by g* —>■ grj' via the anomalous coupling 
of the rj' to two gluons. The calculated rate for B ^ rj' + Xg 
is in rough accord with experiment using a fairly constant glue- 
glue-r/' form factor. This behavior of the form factor may be 
indicative of glueball dominance of the channel. Searches via the 
modes rj'h^h^ {h = tt or K) may be worthwhile. Charmonia 
contributions [i.e B — » rjc,''p{—* rj' + X) + X^ can only account 
for at most 20% of the central value of the signal. Implications 
for B 7]' +Xci and for the corresponding 77 modes are also given. 
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1 Introduction 



Recently CLEO has reported [|T[] a very large branching ratio for the inclusive 
production of rj' (subject to the indicated cut): 

Br(5 + X] 2.2 < Er^, < 2.7GeV) = (7.5 ± 1.5 ± 1.1) x 10^^ (1) 

In this Letter we report on possible mechanisms for this large signal. Two 
interesting origins are: (1) b —>■ s + g* (where g* is an off-shell gluon) followed 
9* ^ v' ~^ 9 due to the anomalous rj'-g-g vertex. (2) b — > charmonia (e.g. 
ipyVc) + followed by ipi^rjc) ^ rj' + X. We discuss these below in turn. 

2 Role of the Anomaly 

We suggest that the dominant contribution to the observed large B 7]' + X 
signal is due to two crucial aspects of the Standard Model (SM). First, the 
SM predicts a very large branching ratio for the charmless 6 — >■ s penguin 
transition p|, ^: 

BriB Xs) ~ 10-2 (2) 

(where Xs is a charmless state containing a s-quark). On the quark level 
these may be regarded as due to 6 ^ sg* transitions, where g* is an off- 
shell gluon, and it is the gluons from this reaction that could be the source 
for the large rj' production. This conversion of gluons to t]' states may be 
accomplished through another important aspect of the SM, namely the QCD 
anomaly: 

d,j^ = ^G^^G^-^ + 2zY: m.q^'q (3) 

q=uas 

where is the singlet axial current]^]. Indeed it is the anomalous contri- 
bution to the divergence of the axial current that distinguishes the SU(3) 
singlet rji state from the octet pseudo-goldstone r^g state. Intuitively it is the 
gluonic cloud of the 7]' which makes it appreciably heavier than the members 
of the octet: vr, K and rjs^. The physical rj and rj' states are related to these 
according to the mixing 

T] = cos Ofjrjs — sin O^rii 

7]' = sin OrjTjs + cos 9j^r]i (4) 
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where the mixing angle ft^ is estimated0 in the range of —10° to —20° so r]' 
is dominantly t]i. 

In order to model the gluonic coupling of the rji we use an ?7'-glue-glue 
effective vertex: 

Hiqlqlql)6''%..pq^q^,e^,e^2 (5) 

where gi, q2 and ei, 62 are the 4-momenta and polarizations of the two gluons 
and a, b are color indices. Thus the picture that we have for the inclusive 
rj' production is a combination of 6 — > sg* followed by g* — >■ rj'g as shown in 
Fig. 1. Indeed there have been many phenomenological attempts to quantify 
the gluonic content of the rj'^, In the effective vertex (^) if is a form 
factor that is a general function of the momenta, qf, gf and q"^,. If we expand 
this function in qf and g| at q^, = m^,, the leading term is generated by 
the QCD anomaly. We will determine H{qf ^ 0, g| ~ 0,m^/) by studying 
ip — > r]''-/. 

As is well known, the b sg* transition can be parameterized in the 
standard model (in the limit 0) by the induced chromo-electric and 

chromo- magnetic form factors 



K7 = +%T.^rs{\V'^)\Fl{e){i,q'-q4)L 



b 



i=uct 



where Vi = V*.Vib, gs is the QCD coupling constant, q = pb — Ps is the 
outflowing gluon momentum, A"^ are the Gell-Mann matrices, R = {1 +7^)/2 
and L = (1 - -f^)/2. 

Using Eq. ^jand Eq. ^we can deduce the differential decay rate for b{pb) — *■ 
v'{qv')^iPs)9{k)- 



dsdt 3847r3m 



F^{2V - {ml - s)W/2) + 2G'^{ml - s)V + AnibFGV 



ij) 

Here s = (g^/ + k)^ = q"^, t = {ps + k)^ and u = {ps + qrj')^, V = {stu — 
ml,mlt)/A, W = -(s -m2,)2/2, G = ^,771^/2/ (2vr2s) and F = fi-nibG. We 
have neglected the small contribution of the w-quark and defined /, = F^—Ff. 
We can now readily calculate the rate for b sgt]' if we assume that 
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fi, H are roughly constant as a function of s = q^: 



r(& - sig) = ^^^^^ifhoix) + |^/./2ri(x) + |^/|t2(x)) (8) 
where Fq = \ Vcb\'^ Gj^ml / {192n^) , x = m^i/ml and Tj(x) is given by: 

ro(x) = xMog(l/x) + — (l-a;)(3-27a;-47a;2 + 13a:^-2x^) 

60 

Tiix) = 2x^(2x + 3)log(l/x) + -(l-x)(l-llx-47x2-3x^) 

6 

r2(x) = -x(3x + 2)log(l/x) + -x)(3 + 47x + llx^ -x^) (9) 

The value of To/Tb where is the total width of the B, may be ascertained, 
at leading order, by factoring the phase space effects of the c-quark into the 
semi-leptonic branching ratio so that if Xc = (mc/mb)^ then for I = e, ii: 

T{B ^ IViX,) = ((1 - 8x, + xl){l - xl) + 12x1 log -)ro (10) 

Xq 

Thus using T{B IViX^)/Tb = 0.1^ we obtain Tq/Tb ~ 0.2. 

Next to deduce H we describe ip — 'yrj' by assuming that the ip is a 
weakly bound state of cc quarks[P] in conjunction with the coupling in Eq. ^ 
[see Fig. 2] . Explicit calculations then show that the amplitude is dominated 
by on-shell gluons leading to 

m ,H cos Ur, — T — \ r-, — 2 U-'-J 

^ 2^2 r(?/> ^ e+e") r^log r 

where r = {irirji/m^Y. 

We thus arrive at if(0,0,m^,) ~ l.SGeV~^ for asijn^) = .25 and 9jj = 
— 17° |1T0|, ^ and other data from 0. As mentioned before, there is consider- 
able uncertainty in the value of 6r]. However, due to our normalization with 
respect to —>■ '-frj', the rate for b —>■ sgrj' is independent of the precise value 
of On- This value of H is in rough accord with the coupling which one expects 



from an anomalous QCD contribution to the g-g-rj' vertex|Q, |TT 



Finally, in order to estimate the rate for b rj'gs we need to know 
the value of Fi and F2. Fi is calculated at leading order, for instance, in 
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where Fi ^ —.20. We can infer that QCD corrections do not effect this result 
greatly by considering next-to-leading order (NLO) calculations of penguin 
T^. To do this we expand in terms of their Oj operators: 



operators m 
C4O4 + CeOe 



(1/2)(C6 - c,)iOe - O4) + (1/6)(C6 + C4)(03 + O5 



+ (C6 + C4) 



(12) 



where Cj are the appropriate Wilson Coefficients of Oi |T^, |13|. We note that 
the last term corresponds to a color octet exchange with the same quantum 
numbers as a gluon. Thus we attribute it to the exchange of a single effective 
gluon from a penguin inferring that Fi = 4:{c4^ + Cq)/ Qs which gives a result 
that agrees with the above to 0(20%). For instance, using the leading order 
calculation in |jl3[ with = 22hMeV so = .21 we obtain Fi = —.168. 
The quantity F2 may be taken directly from this calculation since F2 = Csg = 
.143. 

We can now readily estimate the branching ratio for b — srj'g. Using 
Fi = -.168, Hcoser, = 1.7GeV-^ and Tq/Tb = .2 and taking nib = 4:.8GeV 
we get 1.9 X 10^^ if just the electric form factor is used while the presence of 
the magnetic form factor increases this result by about 50% yielding Br{b — >■ 
sgr]') ~ 2.8 x 10"^ 

The results quoted in [jl| however contain an acceptance cut designed 
to reduce the background from events with charmed mesons[0]. This cut 
requires that the energy of the rj', E^i > 2.2GeV which is equivalent to the 
cut in the recoil mass rrirec ^ 2.35GeV 

In order to understand the rrirec spectrum from a 5-meson, one must 
factor in the fermi motion of the 6-quark with respect to the meson. For this 
purpose, following previous works |l^, we assume that the momentum of the 
b-quark has a probability distribution of the form 



P oc e 



(13) 



which is suggested by a harmonic-oscillator like wave function. Fits to the 
semileptonic decay spectrum at CLE0|]16[ suggest a value of P = 0.287GeV; 
iTif, = 4:.87GeV. In our results we consider P = .15 — AbGeV, though in 
most cases the branching ratios are not strongly dependent on f3. When the 
cut is imposed however, the fraction passing the cut is strongly dependent 
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on mfe. Following CLEO (6 — S7) |T6[ we vary over the range of 4.7 to 
A.dGeV. Using /3 = 0.3 GeV we find that for rrih = 4.7 GeV the branching 
fraction which passes the cut is 6.7 x 10^^ while if = A.dGeV then the 
corresponding branching fraction is 9.8 x 10~^. Furthermore the proportion 
of events passing the cut varies from CcutiE^j') ~ 25 — 33%. 

Thus, taking mi, = A.SGeV and identifying b sgt]' as the dominant 
contribution 111] to the inclusive process we find that after including the ex- 



perimental cuts: 

Br{B T]'X) ^ 8.2 x 10"^ (including cuts) (14) 

which is in rough accord with the experimental result given in Eq. |1|. If we 
include a 150 MeV uncertainty in nib and 0.04 in as then in this model the 
above number varies by about 30%. Fig. 3 shows the rurec distribution for 
P = 0.15, 0.3 and 0A5GeV, = A.SGeV as well as the case /? = O.SGeV 
for rrib = A.GGeV. As can be seen, for a fixed rub, there is some variation 
in these spectra with f3 only in the region rrirec < 1.5G'eV^[0]. The reason 
for the large dependence on nib is readily apparent since as rUb decreases the 
curve shifts to the right and the proportion of events falling below the m^ec 
cut corresponding to E^' = 2.2GeV rapidly drops. 

Although there are several uncertain parameters in the calculation {as, 
rub, /3 . . .) it is quite remarkable that an approximate constant form factor H 
leads to a rate (Eq. (0)) that roughly agrees with experiment (Eq. ([1|)) as 
one of the gluons is considerably off-shell, (ql) ~ 10 GeV^. This behavior of 
the form factor may be an indication that it is dominated by the presence of 
gluonic states in the two gluon channel. The large 77' signal offers a unique 
possibility, in any case, to search for such states through the modes, for 
example, 

Glueball t]' + h+ + (15) 

where h = n or K. Thus it may be worthwhile to closely study the invariant 
mass distribution of (g^/ + qh+ + Qh-Y in the data sample. 

3 Contribution from Charmonia 

Charmonia contributions (See Fig. 4) can be subdivided into three categories: 
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1. The most important of this class of contributions is from b ^ rjc + s 
followed by t]c ^ f]' + X. (See Fig. 4a). The experimental accep- 
tance cut plays a very important role here. It essentially eliminates all 
contributions from a multi-body state X and only a single particle or 
resonance (X) remains viable. Important examples of X are: cr,?7, u, 
V. . . m. We find Ccutimx ~ 2m,) ~ 20%, Ccut{mx ~ m,) ~ 15%, 
Ccuti^x ~ ^uj) ~ 10% and Ccutij^x ~ ^(f>) ~ 1-5%. From the mea- 
sured branching ratios of rjc into such states we estimate that their sum. 



Ei Br{r]c ^ iXi) ~ 10%. Using Br{B ri^ + Xs) ~ 7.0 x 10"^ @, p 
we thus see that, from this mechanism: 



cut 



Br{B ^r]' + X,) = Br{B ^ r/, + X,) x ^ Br{T], r]' + X^) x C 

i 

~ 1.1 X 10-^ (16) 

2. Next we consider b i]c*{^ v') + ^- (See Fig. 4b). We view this 
contribution to arise through a mixing of 7]c and 7]'. Such a mixing is 
predominantly driven by the two gluon intermediate state. It is very 



difficult to see why this mixing angle should be bigger than a degree 
Thus Br{b r]*{-^ t]') + s) = sin^ 6^^ x Br{b ^ r]^ + s) <3x 10"^ 
and is negligible. 

3. Finally there is also the possibility that b ip + s followed hj ifj ^ 
rj' + X. [See Fig. 4a] Again, the acceptance cut effectively requires that 
X be a single particle state. Thus possible examples of X are 7, uj, (j) 
etc. The corresponding branching ratios are expected to be very small, 
each is < 0.5% §, suggesting that EiBriip ^ t]' + Xi) ~ 2%. We 
thus estimate from b ^ ip + s the resulting contribution is Br{B — 
V + X,) ~ 3 X 10"^ 

Thus the main charmonia contribution is via B ^ rjc + Xg, rjc rj' + X 
with an effective Br{B — ^ rj' + Xs) ~ 1.1 x 10"'^ (including cuts), which is at 
most 20% of the observed central value. 

Experimentally, the charmonia contributions coming via cascade decays 
of ip and or rjc can, in principal, be separated by studying the invariant mass 
distributions of the expected final states, such as ri'+rj^r]', (p, u, 7, vrvr, KKetc), 
which should cluster around m^^ and/or m^. 
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Note also that while the rjc — t]' mixing is expected to make a negligible 
contribution it tends to have a distinctive rrirec distribution which is quasi- 
two-body; it peaks around rrirec ~ l-5GeV and diminishes rapidly either side 
of this region. 



4 Other Mechanisms 



In passing we note that there are other mechanisms [|TT]] for making rj' in B 
decays: 

1. There is a 6 — sss penguin process where the s-quark from the gluon 
may pair with one of the s-quarks as shown in Fig. 5a. 

2. There is a 6 ^ suu (or sdd) penguin process where the light quarks 
pair as shown in Fig. 5b. 

3. There is a 6 ^ suu (and also b sdd) penguin process where the u{d) 
can pair with the spectator u{d) to make the r]' as shown in Fig. 5c. 

4. There is a 6 ^ sgg penguin sub-process followed by gg rj' shown in 
Fig. 5d. 

5. The r]' may also be produced through a. b ^ suu tree graph shown in 
Fig. 5e. 



Using the NLO effective Hamiltonian and vacuum saturation we can 
estimate the contributions of Fig. 5a and Fig. 5b. This is straightforward for 
all of the operators except for the contributions of O5 and Og. In these cases 
we need to know 

GrjAli) = V3< 0\q,l^qi\vi > (17) 

Due to the presence of the anomalous contribution there is no direct 
way to deduce the value of Gr^ but a reasonable estimate is perhaps that 
Gfj-^ = Gt, where 

G^ =< 0\u-f^d\n+ >^ mlf\/{mu + m^). (18) 

Using rUu + rrid ~ 15MeV, Fig. 5a together with Fig. 5b give ~ 7.5 x 10~^. 

We estimate the three-body penguin (Fig. 5c) to contribute about 2x 10~^ 
which gives 1 x 10~^ with the cut implemented. In this estimate, we have 
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used as normalization, the CLEO central value for Br{B^ tt I'^vi) = 
(1.8 ±.4 ±.3 ±.2) X 10-^ §. 

Simma and Wyler [|I^] consider gluball production from b —* sgg. Their 
calculation can be adapted to the i]' (Fig. 5d) yielding about 1 x 10~^. 

The b — i> suu tree graph contribution (Fig. 5e) may be estimated to be 
about 5 X 10~® which becomes about 2.6 x 10~^ with cuts. This contribution 
is only to B~^] the B^ decay via the tree is color suppressed. Thus when one 
averages over 5° and B^ this is effectively reduced by a factor of about 2. 

Thus, for b rj'+Xg all of these processes are found to be appreciably less 
than the production through the anomaly coupling (Fig. 1) and indeed not 
even competitive with the charmonia contribution, Eq. (|16D. We estimate 
that these subdominant processes will not alter the main contribution by 
more than about 10%; therefore, for now we can safely ignore them ||1 1|| . 

Inclusive t] production, B —>■ rj + Xg, will in general proceed through the 
processes which we have discussed above for the rj'. The production through 
the gluon mechanism of Fig. 1 is suppressed by approximately tan^(6',,) (up to 
small additional corrections due to the difference in mass). Thus if 6*^ = —17° 
then the total branching ratio with this mechanism is 2.6 x 10~^ while with 
a cut of > 2.2GeV the branching fraction is 8 x 10~^. 

Next we consider the charmonia contribution to rj. From Ref. we see 
that in r]c, ip decays, the branching ratios of final states with t] are roughly 



the same as those with ^7'[§0|- This leads us to suggest that, from charmonia, 
Br{B ^ + X,) ~ 7.3 x 10"^ 

The mechanisms Fig. 5a and 5b for the t] give 1.3 x 10~^. The kinematics 
are such that virtually all of the events will pass the cut. The process shown 
in Fig. 5c gives a total branching fraction of about 3.6 x 10~^ of which about 
1.3 X 10~^ passes the cut. The tree process gives ^ 3 x 10"^ (with cuts) but 
only from B^ decays thus effectively reducing this by a factor of 2 when B^ 
and B^ are averaged. 

It is also interesting to consider the situation when one replaces the 6 — > s 
penguin with a 6 — >■ ci penguin and thus contemplates the inclusive process 
b ^ v'Xd. 

With the exception of the tree graph (Fig. 5e) and the charmonium, all 
of the processes considered above will scale with r^^ = |VJd/V^sp ~ C'(A^) 
where A = sin^c ~ .22 while the charmonium processes will scale like A^. 

The tree graph, on the other hand, leads to T{b duu)/V{b suu) ~ 
1/A^ so the contribution of this process to the branching ratio of 6 — >■ r/'X^ 
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is about 5.4 x 10~^ (including cuts). Therefore this may be expected to 
provide the dominant contribution to r]' + Xa- For 5" this process is 

color suppressed. So from the tree, averaging over the two, B ^ rj' + ~ 
2.7 X 10-5 (with cuts). 



5 Summary and Outlook 

Table |l| summarizes the expected rates, for individual mechanisms, as well 
as the total for the four final states of interest: rj'Xg, Tj'Xd, rjXs and r/Xd- 
For the rj'Xg, the branching ratio is about 3.5 x 10"^ and is about 2.5 x 
10~*, 1 X 10^'^ and 1.4 x 10"^ respectively for the other three, assuming for 
simplicity, = 0.05. The current CLEO acceptance cuts [Eq. (|l|)] 

affect the anomalous contribution (Fig. 1) more severely (about a factor of 
3 to 3.5) than the tree contribution (about a factor of two). As indicated 
in the Table the anomalous graph is the dominant mechanism contributing 
perhaps about 86% of the total for rj'Xg. It contributes about 55%, 44% and 
10% to rj'Xd, rjXs and rjXd respectively. 

An especially notable feature of the numbers in the Table is that the rates 
with cuts for the other modes {rjXd, rj'Xd and r]Xs) are also predicted to be 
quite large, i.e. (0.4 to 1.8) x 10"^, although none is as high as the rj'Xg. 
Their searches clFG Gel gerly awaited as comparison with experiment will be 
very instructive. 

In closing, it is important to note that the h — > srj'g mechanism should be 
regarded as the dominant contributor to the multiparticle inclusive process 
B —>■ rj' + Xg, where Xg is a kaon together with at least one pion. It is 
expected [0 to make a negligible contribution to the exclusive two body 
mode {B ^ 7]' + K) . We recall that experimentally B ^ rj' + K seen [jl| to 
be about 10% of the inclusive signal B ^ t]' + Xg. To facilitate contact with 
theoretical models, it may be better if the two body signal is not included in 
the inclusive rate, Eq. |[ 

It is extremely interesting that the rj' becomes such an effective probe of 
the b — > sg* process. It has long been known that these charmless final states 
have a very large branching ratio @]. The rj' could therefore become a very 
powerful tool for gluonia searches and for illuminating our understanding of 
the penguin dynamics, and, in particular for the quest for direct CP violation 



as the same final states are expected to be rich in CP asymmetries [pi 
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Figure Captions 



Figure 1: The main mechanism for the inclusive 77' production: b — > sg* 
graph followed by g* — > grj' where the black circle denotes the anomalous 
ggrj' couphng. 

Figure 2: — > ^ygg followed by gg — > rj' via the anomalous gg — > rj' 
vertex. 

Figure 3: Recoil mass distribution for B Xsi]'s from the sub-process 
b — > gTj's. The events falling to the left of the vertical line pass the -E^' > 
2.2GeV cut. For = 4.8 GeV, the solid curve shows the distribution 
for P — O.SGeV while in the region m^ec < 1.6GeV the distributions for 
P = O.lbGeV (dashed) and (3 = QAbGeV (dotted) are also shown. For 
mb = 4.6 GeV and f3 = 0.3 GeV the rrirec distribution is shown by the dot- 
dashed curve. Each of the above curves is normalized so that its integral is 
unity. 

Figure 4: Charmonia contributions: a) b ^ ip,ric + s followed by ip, rjc —>■ 
Tj' + X. h) b ^ T]') + s; X denotes rjc rj' mixing. 

Figure 5: Feynman diagrams for the sub-dominant contributions to 5 — >■ 
Tj' + Xg. The four quark penguin operators (i.e. O3 — Oe) are represented 
by hexagons. The diagrams shown are: (a) the b — > sss process where 
the s-quarks may combine with either of the s-quarks to form an r)'; (b) the 
b — > suu or sdd process where the uu (dd) forms the 1]'; (c) A typical penguin 
graph producing two gluons that fuse to form a rj'; (d) a 6 — suu or sdd 
process where the light quark combines with the spectator; (e) the b — > suu 
tree process. 
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Table 1: Contributions from different meclianisms to tlie Br(B — >■ f]'{i]) + 
Xg^d) in units of 10^^. The number from eacli graph is given along with the 
total, without cuts. The effect of CLEO acceptance cuts on the total is also 
indicated. Note = 20\Vtd/Vts\^. 



Mechanism 


Type 




v'Xd 


vXs 


vXd 


Fig. 1 


Anomaly 


280 


lArl 


26 


l.2,rl 


Fig. 4 


Charmonia 


70 


3.4 


70 


3.4 


Fig. (5a + 5b) 


penguin-4q-op 


0.75 


0.04r2 


1.3 


0.07r^ 


Fig. 5c 


penguin-4q-op 


2.1 


O.lr^ 


3.6 


0.18r2 


Fig. 5d 


penguin-2g-op 


1.0 


0.05r^ 


0.1 


0.005r2 


Fig. 5e 


tree 


0.25 


5.2 


0.4 


8.6 


Total (no cuts) 




354 


14r2 + 8.6 


101 


1.6r2 + 12 


Tc)tal {with cuts) 




96 


1.2/-^ + 3.2 


17.6 


.5 1/-5 + 3.3 


Anomalous ^^j^j^ ^^^g^ 




~ 86% 


~ 55%"^ 


~ 44% 


~ io%«) 



a) assuming r^ — l. 
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Figure 1 




Figure 2 
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Figure 4 




(4b) 



Figure 5 




